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We demonstrate that a unidirectional coupling between a pattern forming system and its replica 
induces complete synchronization of the slave to the master system onto a spatiotemporal chaotic 
state. 



In recent years, synchronization of complex systems 
have attracted a great interest in the scientific commu- 
nity [lL as wei l as m the literature oriented to lay audi- 
ences y. This indicates the behavior of two (or many) 
systems (either equivalent or non equivalent) that adjust 
a common feature of their complex dynamics due to a 
coupling or to a forcing. 

For time chaotic systems, four types of synchroniza- 
tion have mostly been studied, namely complete syn- 
chronization (CS), phase (PS) and lag (LS) synchroniza- 
tion, and generalized synchronization (GS). CS refers to 
a process whereby two interacting systems perfectly link 
their chaotic trajectories, thus remaining in step with 
each other in the course of the time Jj. GS implies the 
hooking of the output of one system to a given function 
of the output of the other system Q. PS is character- 
ized by a locking of the phases of the two signals, also 
in the absence of a substantial correlation between the 
two chaotic amplitudes Finally, LS consists in the 
hooking of one system to the lagged output of the other 
@. All these effects have been explored in natural phe- 
nomena 0, and laboratory experiments and unified 
approaches to describe and measure synchronization 
states have been proposed. 

When the interest shifted to space-extended systems, 
synchronization phenomena were shown in large popula- 
tions of coupled chaotic units and neural networks |lfj| . 
globally or locally coupled map lattices [llj , and pattern 
forming systems governed by partial differential equa- 
tions [12| . Here, however, all theoretical and numeri- 
cal progresses were accompanied by a substantial lack of 
experimental verifications. Precisely, CS of spatiotem- 
poral patterns was first observed in chemistry for 
two mutually coupled Belouzov-Zhabotinski cells, where, 
however, the resulting synchronized state corresponded 
to the suppression of spatiotemporal complexity and the 
emergence of a common spiral behavior. Later on, LS 
was observed in a pair of unidirectionally coupled non- 
linear optical systems |l4j . The evidence here was given 
in terms of an improvement in the lagged correlation be- 
tween the master and slave patterns. Finally, GS was 
demonstrated in an open loop liquid crystal light modu- 
lator with optoelectronic feedback Il5j by the use of the 
so called auxiliary system method |4j. 

In this Letter we report the first direct experimental 
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FIG. 1: a) Experimental setup. Main loop: an expanded 
laser beam is closed through a non-linear Kerr-like medium 
(liquid crystal optical valve). BS: beam splitters; A: aperture 
in a Fourier plane. Control arm: PC: personal computer; 
LCD: liquid crystal display. 



evidence of complete synchronization on unidirectionally 
coupled pattern forming systems. At variance with what 
observed in Ref. [I3> the resulting synchronized state 
here corresponds to a spatiotemporal chaotic dynamics 
where the slave system is identically attained to the be- 
havior of the master system. 

The experimental setup is sketched in Fig. 1. A main 
optical feedback loop (MOFL) consists of a Liquid Crys- 
tal Light Valve (LCLV) with optical feedback 16]. The 
LCLV induces on the reading light a phase delay pro- 
portional to the writing intensity (Kerr-like effect). This 
proportionality relation holds for all experimental param- 
eters used in the present investigation. The pattern form- 
ing mechanism acts on an homogeneous pump wave of in- 
tensity I p sent onto the front face of the LCLV. The wave 
is totally reflected and acquires a spatial phase modula- 
tion. Diffractive propagation along the MOFL provides 
conversion of phase into amplitude modulations. Due 
to the Kerr-like behaviour of the LCLV, these amplitude 
modulations are converted back to phase modulations, so 
that eventually a positive feedback establishes for some 
spatial frequencies, which are destabilized resulting in 
pattern formation. 

An additional electro-optic control loop is constituted 
by a video-camera and a personal computer driving a liq- 
uid crystal display (LCD). The control signal is a laser 
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beam that traverses the LCD before being injected into 
the MOFL. The LCD display is operating in transmis- 
sion, and encodes linearly the gray level images output 
by the PC, onto the laser beam traversing it. 

When the control loop is open, the dynamics of the 
optical beam phase 6(x, y, t) at the LCLV output can be 
described by [H[l| % = -±(0 - <fo) + DV 2 (j) + al fb , 
where </>o is the working reference phase, r the LCLV 
relaxation time, D a diffusion coefficient, a the LCLV 
nonlinearity strength, and Ifb(x, y, t) the feedback inten- 
sity at the input plane of the FB. In is a nonlinear (and 
nonlocal) function of the phase <j> [ijl UJ1 • 

With increasing the pump intensity I p above the 
threshold ithr of pattern formation, the homogeneous so- 
lution destabilizes, and an hexagonal pattern arises. This 
allows to introduce a reduced pump parameter I = j 2 - . 
A further increase in / above unity leads eventually to a 
destabilization of hexa gons in favor of a regime of space- 
time chaos (STC) |l7l ll9j- An aperture A (located in 
a Fourier plane of the MOFL) has the role of limiting 
the spatial frequency bandwidth B of the system, which 
is another control parameter. Troughout the experiment 
here reported, B is kept fixed at 1.5 times the diffractive 
spatial frequency of the system. This is the frequency of 
the hexagonal pattern which bifurcates at the threshold 
for pattern formation. 

When the control loop is closed, a fraction of Ifb is ex- 
tracted and detected by a video-camera, which is inter- 
faced to a personal computer (PC) via a frame grabber. 
The PC processes the input image, and sends a driving 
signal to the LCD, upon which a plane beam of intensity 
io incides. The transfer function T(x,y,t) of the LCD is 
the sum of a constant mean transfer coefficient To and a 
modulation signal s(x,y,t), which we set to be propor- 
tional to the error signal between the actual pattern in- 
tensity I ft, and a desired time dependent target pattern 
Ir{x,y,t) [s(x,y,t) = -(7/I0) (If b (x,y,t) - I T (x,y,t))]. 
Further real time processing performed by the PC in- 
cludes the evaluation of s, and the calculation of the 
cross-correlation between I ft, and It- The resulting ac- 
tualization time for s is of the order of 200 to 300 ms, 
to be compared with the characteristic time of the pat- 
tern dynamics (of the order of 1 sec. for the parameters 
used in our experiment). The diffractive scale of the sys- 
tem is V2AL ~ 300 fim (A = 514 nm being the laser 
wavelength, and L = 90 mm the free propagation length 
in the MOFL). On the other side, the control area is 
~ 2000 x 2000 ^m 2 , and the control signal is made of 
128 x 128 pixels. This grants us a spatial resolution of 
~ 20 pixels per typical pattern wavelength. 

With the help of such real-space real-time control tech- 
nique we have recently given evidence that two dimen- 
sional stationary target patterns with arbitrary symme- 
tries and sha pes can be effectively and robustly stabilized 
within STC |20j. Here, instead, we aim to demonstrat- 
ing complete synchronization in a unidirectionally cou- 
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FIG. 2: a) Snapshot of the master dynamics recorded with 
open control loop at I = 3. a) Space (vertical)-time (hori- 
zontal) dynamical evolution of a line of (a) the output free 
running signal from the MOFL (master dynamics); (b) the 
slave dynamics in conditions of closed control loop; (c) the 
synchronization error, corresponding to the difference of (a) 
and (b). The gray scale in (c) is such that black corresponds 
to negative values, white to positive values and gray to values 
close to 0. 

pled scheme between two identical systems in a regime 
of STC. For this purpose, we initially let the control loop 
open and record (over a time interval T) the free evolu- 
tion of the system for a value of / at which the uncon- 
trolled dynamics displays STC. A qualitative character- 
ization of the resulting dynamics has been given in |20| . 
The signal consists basically of a set of closely packed 
diffractive spots, evolving in space and time in an unpre- 
dictable way. 

After the registration of this dynamics, which we re- 
fer to as the Master Dynamics (MD) henceforth, further 
free evolution is granted to the uncontrolled system, so 
that after a few seconds the configuration of the MOFL 
output is totally uncorrelated with the initial frame of 
the MD. At this point, we close the control loop and re- 
play the registered MD as the target pattern Ix{x,y,t). 
In this way, we are implementing a unidirectional cou- 
pling scheme between two identical systems starting from 
fully uncorrelated initial conditions. By repeating the re- 
playing procedure of the MD with increasing values of 7 
(hereinafter called the coupling strength parameter), we 
eventually observe full synchronization between the con- 
trolled output of the MOFL [the slave dynamics (SD)] 
and the MD. 

Complete synchronization between SD and MD is 
shown in Fig. 2, reporting the space (vertical)-time (hor- 
izontal) dynamical evolution of the central vertical line of 
pixels for the master dynamics (2a), slave dynamics (2b) 
and the difference between the two, (2c). Fig, 2 is taken 
at 7 = 0.8, and shows how the SD closely follows the MD 
at any time during control. As visible in fig, 2c, the syn- 
chronization error E(x,y,t) = Ifb(x 7 y,t) — lT(x,y,t) is 
close to vanishing nearly always and everywhere in time 
as a result of the complete synchronization process. 

Notice that the final CS state is here realized within a 
full STC regime, at variance with what reported in Ref. 
13] for a bidirectional coupling between two excitable 
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FIG. 3: Experimental point distributions in the (isd,imd) 
plane (see text for definition of the variables). The syn- 
chronization manifold is represented by the diagonal line 
£sd = imd. (a) I = 2, 7 = 0.8, b) I = 3, 7 = 0.8, c) 
1 = 2, 7 = 0, d) 7 = 3, 7 = 0. 



media, where the emergence of synchronization was as- 
sociated with the suppression of space-time chaos in the 
system. 

An independent way of visualizing the emer- 
gence of CS is to pick randomly a set of N 
points [(xi,yi); i = 1,...,N] in both the MD 
and the SD, and to plot the variable Xsr>(t) = 
{Ifb(xi, Ui,t), i = 1, . . . , N} vs. the corresponding vari- 
able ccmd(^) = {Ir{xi, yi, t), i=l,...,N}. The more 
the distribution of points in the (xsd, imd) plane ap- 
proaches the synchronization manifold (the diagonal line 
xsd = imd)i the cleaner a CS state is set in our system 
at all times. Experimental results are reported in Fig. 3 
for two values of the coupling strength 7. Precisely, Fig. 
3a (3c) shows the distribution of points in the (xsd, xmd) 
plane for I = 2 and 7 = 0.8 (7 = 0), while Fig. 3b (3d) 
refers to the same situation for I = 3 and 7 = 0.8 (7 = 0). 
In both cases, it is apparent that increasing the coupling 
strength induces a point distribution much closer to the 
diagonal line than the uncoupled case, giving evidence 
that a CS state has arisen in the experiment. 

A more quantitative measurement of CS can be given 
by monitoring the behavior of the time dependent cross- 
correlation function C(t) =< i/&(r,t) ■ 7r(r,t) >r be- 
tween the instantaneous patterns in the SD and MD dur- 
ing the synchronization process (< . . . >r denotes here a 
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FIG. 4: Cross correlation function C(t) (see text for defini- 
tion) vs. time during some synchronization trials for 1 — 3. 
a) 7 = 0.24, b) 7 = 0.40, c) 7 = 0.8. In all cases t = 
indicates the instant at which the MD starts to be replayed 
in the control loop. 



spatial average in the plane r = (x,y)). C(t) is by defini- 
tion vanishing for linearly uncorrelated systems, whereas 
C{t) ~ 1 for fully synchronized dynamics. 

Fig. 4 reports the temporal behavior of C(t) for 1 = 3 
and for three different values of the coupling strength [a) 
7 = 0.24, b) 7 = 0.4, and c) 7 = 0.8]. In all horizontal 
axes, t = has been taken as the instant at which the 
MD starts to be replayed in the control loop. A first im- 
portant observation is the cross correlation starts from 
a non vanishing value at t = 0. This is because the 
uncontrolled STC dynamics has a non zero mean field, 
as it can be appreciated from inspection of Fig. 2a-2b. 
Namely, the uncoupled MD and SD have a certain degree 
of "phase rigidity", i.e., even if there are chaotic fluctu- 
ations, bright (dark) areas remain more or less bright 
(dark) for most of time. Similar properties have been 
observed experimentally and discussed in various other 
cases of space extended systems giving rise to STC dy- 
namics [2lj |. 

At low values of the coupling strength (7 = 0.24), 
no synchronization is set in the system. This is visi- 
ble in Fig. 4a, where C(t) experiences large fluctuations 
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in time around a mean value not substantially different 
from the initial correlation level. For intermediate cou- 
pling strengths (7 = 0.4 in Fig. 4b), a partial synchro- 
nization emerges in the system after a transient time, 
though several deviations of the SD from the MD still 
remain, reflected by the rather large fluctuations around 
the asymptotic value of C(t) visible in Fig. 4b. Finally, 
a further increase in 7 induce a full CS in the system 
(see the case 7 = 0.8 in Fig. 4c), where the asymptotic 
value of the cross correlation approaches unity, and the 
residual fluctuations in C(t) shrink considerably. Notice 
that, as 7 increases, the transient time before reaching 
CS decreases. The global picture depicted in Fig. 4 con- 
firms that CS here is a threshold phenomenon, as it was 
introduced originally for time-chaotic systems Q ■ 




FIG. 5: Time average cross correlation function C (see text 
for definition) vs. the coupling strength 7 for 1 = 2 (circles), 
1 = 3 (squares), and 1 = 4 (triangles). In all cases, the solid 
lines connecting points are drawn as guides for the eyes. 

The global scenario of observed CS is illustrated in Fig. 
5, where the time average of the cross correlation function 
C =< C(t) >t (< . . . >t here denotes a further time 
average over the full time interval T) is reported vs. the 
coupling strength at different values of pump intensities 
/. In all cases, our proposed coupling scheme is effec- 
tive in inducing CS, as C considerably increases with 7 
with respect to the corresponding uncoupled values. An 
increase in the pump intensity / leads to a progressive de- 
terioration of CS, reflected by smaller asymptotic values 
of C{t). 

In conclusion, we have demonstrated that a unidirec- 
tional coupling between two identical pattern forming 
systems induces complete synchronization of the slave to 
the master dynamics, with a resulting synchronized be- 
havior corresponding to a spatiotemporal chaotic state. 
Though realized with a nonlinear optical experiment, the 
coupling scheme used (based on a real-space real-time 
control of a recorded and replayed space-time chaotic dy- 
namics) can in principle be implemented in many other 



physical and chemical pattern forming systems. Further- 
more, the great flexibility offered by the proposed cou- 
pling scheme can be exploited to drive the slave system 
onto a generic desired dynamics with arbitrary symme- 
tries and shapes in space, as well as arbitrary behaviors 
in time. 

Work partly supported by MIUR-FIRB project n. 
RBNE01CW3M-001. 
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